N = [2:65,70:5:200]';

x = load('tableData.dat');

Hnumerical = x(:,2) ./ 2;

o = optimset('MaxFunEvals',4000,'MaxIter',4000,'Display','off');

N = [2:65,70:5:200]';

x1200a = lsqcurvefit(@fitFunctionForH,[0,0,0,0,0,0],N,Hnumerical,[],[],o);
x2200a = lsqcurvefit(@fitFunctionForHApprox,[0,0,0,0,0,0],N,Hnumerical,[],[],o);
x1200b = lsqcurvefit(@fitFunctionForH1415,[0,0,0,0,0,0],N,Hnumerical,[],[],o);
x2200b = lsqcurvefit(@fitFunctionForHApprox1415,[0,0,0,0,0,0],N,Hnumerical,[],[],o);


N = [2:65]';
n = [2:65,70:5:120]';

x165a = lsqcurvefit(@fitFunctionForH,[0,0,0,0,0,0],N,Hnumerical(1:64),[],[],o);
x265a = lsqcurvefit(@fitFunctionForHApprox,[0,0,0,0,0,0],N,Hnumerical(1:64),[],[],o);
x165b = lsqcurvefit(@fitFunctionForH1415,[0,0,0,0,0,0],n,Hnumerical(1:length(n)),[],[],o);
x265b = lsqcurvefit(@fitFunctionForHApprox1415,[0,0,0,0,0,0],N,Hnumerical(1:64),[],[],o);

fprintf('\r\n \r\n Using (1-log2) For Approx Fit, 65 then 200 \r\n');
fprintf('%f ',x265a);
fprintf('\r\n');
fprintf('%f ',x2200a);
fprintf('\r\n');

fprintf('\r\n \r\n Using (1-log2) For Exact Fit, 65 then 200 \r\n');
fprintf('%f ',x165a);
fprintf('\r\n');
fprintf('%f ',x1200a);
fprintf('\r\n');

fprintf('\r\n \r\n Using (14/15-log2) For Approx Fit, 65 then 200 \r\n');
fprintf('%f ',x265b);
fprintf('\r\n');
fprintf('%f ',x2200b);
fprintf('\r\n');

fprintf('\r\n \r\n Using (14/15-log2) For Exact Fit, 65 then 200 \r\n');
fprintf('%f ',x165b);
fprintf('\r\n');
fprintf('%f ',x1200b);
fprintf('\r\n');
